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Abstract 

Large assemblies of nonlinear dynamical units driven by a long-wave fluc¬ 
tuating external field are found to generate strong turbulence with scaling 
properties. This type of turbulence is so robust that it persists over a finite 
parameter range with parameter-dependent exponents of singularity, and is 
insensitive to the specific nature of the dynamical units involved. Whether or 
not the units are coupled with their neighborhood is also unimportant. It is 
discovered numerically that the derivative of the field exhibits strong spatial 

intermittency with multifractal structure. 
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We report our discovery of a new type of turbulent behavior which arises generally 
in large assemblies of simple dynamical units driven by a long-wave randomly fluctuating 
held. The driving held may actually be a self-generated internal held due to long-range 
interaction, and this particular situation was studied in a previous paper [1] where a rough 
explanation of the origin of power-law correlations was also given. The present study thus 
aims at expanding as much as possible the class of systems capable of exhibiting the same 
type of turbulence, and also proposing a more transparent and coherent explanation of the 
phenomena. 

An illustrative example is given by an array of uncoupled logistic maps f(X) = aX(l— A") 
with driving 

X n +i(j) = f(X n (j )) + h n (j), j = (1) 

where h n (j ) = y (l + cos[ 27 t{-^ + t/’n}]), VVi being a random variable in the interval (0,1] 
with uniform distribution. Were it not for spatial dependence of h n , Eq. (1) would represent 
N identical copies of a randomly driven map; the dynamics of such an ensemble was studied 
in [2], Making h n nonuniform changes the problem completely. Let the parameter values 
be set such that individual maps are entrained to h n in the sense that their maximum 
Lyapunov exponent (common to all maps) is negative. LInlike h n , however, the corresponding 
amplitude profile is not smooth at all, a typical example of which is displayed in Fig. la. 
Such ill-behaved nature of the pattern is even amplified in Fig. lb which shows a strongly 
intermittent pattern of the differential amplitudes Y(j ) = \X(j + 1) — X(j)\ constructed 
from Fig. la. 

Before proceeding to further numerical study, some theoretical predictions will be made 
as to the statistics of turbulence to be shared by the above system or more general assemblies 
of units under long-wave random driving. For this purpose, it seems more convenient to work 
with a picture in which the dynamical units form a quasi-continuum rather than a lattice, 
and the driving held has a characteristic wavelength of 0(1). Our primary concern is to 
understand how a simple driving held alone can generate a nontrivial correlation between 
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uncoupled units. 

Suppose that each unit is represented by a discrete-time dynamical system X n+1 = 
f(Xn). Let the units be driven by an additive random force h n which is smooth in space 
and statistically invariant with respect to spatial translations. The unit at site r is governed 
by the equation 

X n+ i(r) = f(X n (r)) + h n (r). (2) 

Analogously to fully developed fluid turbulence [3], let us consider various moments of 
the amplitude difference between two sites. We thus concentrate on a pair of units at 
sites r 0 and r 0 + x with distance x = \x\ satisfying x <C 1. The amplitude difference 
y n {x) = X n (r 0 + x) — X n (r 0 ) obeys the equation 

y n +1 = L n y n + 0{\y n \ 2 ) + A h n , (3) 

where L n y n is the linearization of f(X ) about X = X n (r 0 ), and A h n = h n (r 0 +x) — h n (r 0 ) 
is a quantity of 0(x). Equation (3) describes a multiplicative stochastic process [4] with small 
additive noise. Similar equations have recently aroused considerable interest in connection 
with on-off intermittency and related phenomena [5]. While on-off intermittency refers to 
a certain type of temporal self-similarity peculiar to a special parameter value, our major 
concern below is a spatial self-similarity observable over an open parameter range. 

Equation (3) may be simplified by neglecting all eigenmodes of L n other than the least 
stable one. This leads to a scalar equation for y n = \y n \ 

y n + 1 = e Xn y n + 0(yl ) + b n x , (4) 

where X n is the local Lyapunov exponent of the unit at site r 0 , and b n is a randomly changing 
factor of 0(1). For sufficiently small x , there is a range of y satisfying x -C y -C 1 where 
both the nonlinear and inhomogeneous terms in (4) are negligible. We are thus left with a 
linear equation y n +i = e Xn y n or z n+ \ — z n = X n in terms of a new variable z n = In y n . If the 
random process of X n is Markoffian, which we assume, the probability density Q n (z) for z n 
evolves in this linear regime according to 
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( 5 ) 


Qn-\- 1(^) / w(X)Qn(z A)rfA, 

J—oo 

where w(A) is the normalized probability density for A n . Equation (5) admits a stationary 
solution of the form Q(z) oc exp (/3z) = 'if . Thus, the corresponding probability density for 
y n , denoted by P(y), becomes 

P(y ) oc y~ l+ f (6) 

where f3 is determined as a nontrivial (i.e., nonzero) solution of 

/ OO 

e~ f3x w(X)dX = 1. (7) 

-OO 

Note that for sufficiently small /3, we have 

P = 2A/A 2 , (8) 

where the bar means the average with respect to w( A). We have now to modify (6) by taking 
into account the effects of the nonlinear and inhomogeneous terms in (4). The nonlinearity, 
which is assumed to work in such a way that the unstable growth of y n be saturated when 
A„ > 0, may roughly be incorporated by introducing a cutoff in P(y) at y = 1. On the 
other hand, the inhomogeneous term will come into play when y n becomes 0(x ) or smaller, 
thus suppressing the power-law divergence of P(y) there. For the purpose of qualitative 
argument, one may therefore use the following simple model for P(y): 

P(y ) = Cx~ 1+f3 (y < x ), Cy~ 1+l3 (x < y < 1), 

0 (y > 1), C : normalization const. (9) 

This form allows us to calculate the q-th moment < y(x) q > for arbitrary q. For simplicity, 
only positive values of q will be considered below. In the subcritical regime (/? < 0, i. e., 
A < 0), where the dynamical units are entrained to the driving field, we obtain power-law 
moments 

< y(x) q >~ x q (q < |/3|), x m (q > \f3\). (10) 
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The result of q- independent exponent valid for higher moments (q > |/3|) is anomalous 
reflecting strong non-Gaussianity of P(y). In the postcritical regime (/? > 0), all moments 
possess an x-independent part, while the residual part still obeys a power law: 

< y(x) q >~ P(q + /3) _1 + 0(xP). (11) 

The reason why the amplitude difference in the postcritical regime is nonvanishing for van¬ 
ishing x is that the two units in question have lost their respective synchrony with h n , 
implying also the loss of their mutual synchronization. Note that (10) and (11) are asymp¬ 
totic formulae valid for x — > 0 under fixed p. Near \P\ = q and 0 under fixed x, however, 
there exist crossover regimes (Cl) | (P + q) hi x <C 1 and (C2) |/31nx| <C 1, respectively, in 
each of which we have < y(x) q >~ x^\ lnx| and | lnx| -1 . 

A few more remarks are now given on the cases of q = 2 and 1 for which our theory 
recovers our previous results [2], We obtain from (10) and (11) the second moment 

< y(x) 2 > ~ x 2 (P < — 2), x^ (—2 < P < 0), 

m +P)- 1 + 0(x 11 ) (13 > 0), (12) 

while in the aforementioned crossover regimes, we have < y(x) 2 >~ x 2 |lnx| (Cl) and 
1/| lnx| (C2). 

The case q — 1 is related to the length of an amplitude versus space curve. This is 
because the length S(x) for the part of an amplitude profile contained in the unit interval, 
when measured with the resolution of the minimum length scale x, is given by S(x) ~ 
x _1 < y(x) >. Applying (10) and (11), we thus obtain S(x) ~ const. (/? < —1), x^~ l 
(—1 < P < 0), and x^ 1 (P > 0). In the crossover regimes, however, these must be replaced by 
S(x) | In x\ (Cl) and l/(x| lnx|) (C2). The fractal dimension Df defined by S(x) ~ x l Df 
thus becomes 

D f = 1 (P<- 1), 2-|/3| (—1 </3<0), 2 (/3>0) (13) 

except for the crossover regimes. 

The above arguments on discrete-time dynamics can easily be carried over to continuous- 
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time dynamics. One needs only make replacements n —> t, n + 1 —> t + dt and A n —> A (t)dt. 
Then, (4) becomes y = A (t)y + 0(y 2 ) + b(t)x, and (5) reduces to a Fokker-Planck equation 
Q = —A d z Q + \\ 2 d 2 z Q. The latter admits a stationary solution Q(z) oc exp (fiz), and the 
corresponding P(y) is the same form as (6) with f3 given by (8). 

In order to test the validity of our argument, the array of logistic maps (1) has been 
analyzed numerically. In Fig. 2, we display P(y) versus y for some values of K , with x fixed 
at a sufficiently small value. As expected, P(y) exhibits a power-law dependence on y for 
not too small or too large y , with the exponent depending on K . Figure 3 shows moments 
< y(x) q > versus x for some values of q. Their power-law dependence on x is clear, but 
the observed change of the exponent with q, indicated in the small box, is not so sudden 
across q — \(3\ as the theory predicts. As a possible source of this discrepancy, except for 
the existence of the crossover regime Cl, the assumed form of P(y) in (9) might be too 
simplistic, especially in introducing a sharp cutoff at y — 1. 

Up to this point, we have considered uncoupled units. We now show some evidence that 
inclusion of short-range coupling leaves the above-described power-law behavior of moments 
essentially unchanged. As an illustration, we modify (1) with additional diffusive coupling: 

x„ + 1 (j) = f(x n (j)) + Mi) 

+§{f(X„(j + 1)) + f(X n (j - 1)) - 2 f(X n (j))j. (14) 

Without the forcing term h n , the above model would be identical with the usual coupled 
map lattice. In Fig. 4, the second moments < y(x) 2 > are compared between the two 
systems, one with diffusive coupling (D = 0.1) and the other without. The deviation from 
a power law in the presence of coupling is limited to the range covering ten or so units out 
of N (= 4096). This defines a lower cutoff length Xd similar to the dissipation length in 
fully developed fluid turbulence. Although Xd will increase with D like Xd oc \/~D , we have 
a prefactor N _1 , so that Xd can be made arbitrarily smaller than 1 (i.e., the upper cutoff) 
by increasing N indefinitely. Thus, the intermediate range of x , which is similar to the 
inertial subrange, has a sufficient extention over which the power-law nature of correlations 
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is practically unaffected. 

It is also worth noting that whether the dynamical units involved are themselves chaotic 
or not is unimportant to the power-law nature of turbulent fluctuations. This is illustrated 
by driven rotators 

4>j — 1 — ccos <pj + hj(t), |c| < 1. (15) 

Obviously, each unit by itself could never be chaotic. The driving field is assumed to be 
of the form hj(t ) = K cos[27t{^ + ip(t)}], where we let ip(t) behave like the position of a 
Brownian particle, i.e., ip = av and v = —v + £(t) with £(t) representing Gaussian white 
noise. Our assertion is confirmed by Fig. 5 which demonstrates power-law dependence of 
the second moment < y(x) 2 > on x with parameter-dependent exponent. 

Some new aspects of our turbulent held are revealed through an analysis of the differential 
amplitudes Y(j) or quantities defined similarly when the spatial dimension is two or higher. 
We call such a held the F-held. The situation is analogous to fully developed fluid turbulence 
where the study of the energy dissipation held provides rich information which would hardly 
be available from the study of the velocity held alone. We will restrict our discussion to ID 
systems below. Note that in the absence of short-range coupling a true derivative dX/dx 
may not exist in the continuum limit, especially when the amplitude profile is fractal. The 
Y -held must then be redefined as Y(x) = <5 -1 |A^(a; + 5) — X(x)\ with finite but sufficiently 
small 5. 

Spatial intermittency of Y{x) as exemplified in Fig. lb may be analyzed similarly to the 
case of on-off intermittency. This is achieved by measuring the probability density p(l) for 
the space-interval l over which the units are in the laminar state, namely, their Y -values 
stay below a certain threshold Y 0 . Such an analysis was done for the driven logistic maps 
(1). It is clear from Fig. 6 that, as in the on-off intermittency, p{l) exhibits an inverse power 
law. We confirmed that, for not too small or too large Yq, the exponent is insensitive to the 
choice of Y 0 , but depends on K. 

A more thorough characterization of the Y -held is provided by the generalized fractal 
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dimensions D q = (q — 1) _1 lim e ^ 0 In J2i lAI l n e [6]- Here, the measure /q of the i-th box of 
size e is dehned as being proportional to the integral of Y(x) within the same box, with the 
condition of the total measure being normalized. Similar multifractal analysis was performed 
for the energy dissipation field of fully developed fluid turbulence [7]. Figure 7 shows D q 
obtained for the driven logistic maps (1). Note that D 0 = 1, which is simply because Y(x) 
is nonvanishing almost everywhere. An open problem is to explain the various singularities 
of the Y -held and relate them to the singularities of the amplitude field. 

In conclusion, the type of turbulence discussed in this report is so robust that similar 
phenomena should exist quite universally. They may appear in a wide variety of coupled 
and uncoupled systems such as reaction-diffusion systems, fluids and biological populations, 
once placed in a long-wave randomly fluctuating environment. 

The authors thank P. Marcq for fruitful discussions and careful reading of the manuscript. 
The present work has been supported by the Japanese Grant-in-Aid for Science Research 
Fund from the Ministry of Education, Science and Culture (No. 07243106). 
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FIGURES 

FIG. 1. (a) Instantaneous amplitude profile for the driven logistic maps (1). N = 1024, 
K = 0.2, a = 3.7(1 — K). (b) Profile of differential amplitudes Y(j ) constructed from (a). 

FIG. 2. Probability density P(y ) in logarithmic scales for the driven logistic maps (1). For 
each K value, P(y ) obeys a power-law in the intermediate range of y. The exponent changes as 
-1.40, -1.77 and -2.33 with increasing K. a = 3.7(1 — K ), x = 1024^ 1 . 

FIG. 3. Moments < y(x) q > vs. x in logarithmic scales for the driven logistic maps (1), showing 
power-law dependence on x for each q with ^-dependent exponent i](q). Numerical r/ vs. q curve is 
displayed in the small box. K = 0.2, a = 3.7(1 — K). 

FIG. 4. The second moment < y(x) 2 > vs. x in logarithmic scales for the driven logistic maps 
(1) with additional diffusive coupling ( D = 0.1) and without ( D = 0). K = 0.2 a = 3.7(1 — K ), 
N = 4096. 

FIG. 5. The second moments < y(x) 2 > vs. x in logarithmic scales for the driven rotators 
(15) with different values of c. In numerically integrating (15), a simple Euler scheme is adopted 
where f(t ) is replaced with uniform random numbers over the interval (—At, At] given at each 
fundamental timestep of At = 0.05. K = 1.0, a = 30.0. 

FIG. 6. Size distribution p(l ) of the laminar domains of the 

T-field (= Y(x) = \X(x + 5) — X(x)\/5) in logarithmic scales for the driven logistic maps (1). 
K = 0.2, a = 3.7(1 - K), N = 1024, 5 = A -1 , Y 0 = 0.2, 0.4. 

FIG. 7. Dimension spectrum D q of the Y -field for the driven logistic maps (1). K = 0.2, 
a = 3.7(1 - K), N = 1024, 5 = N~K 
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